A Modified Generalized Chaplygin Gas as the 
Unified Dark Matter-Dark Energy Revisited 



Xue-Mei Deng 

Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China 
E-mail: xmd@pmo .ac.cn 

Abstract. A modified generalized Chaplygin gas (MGCG) is considered as the 
unified dark matter-dark energy revisited. The character of MGCG is endued with 
the dual role, which behaves as matter at early times and as an quiessence dark energy 
at late times. The equation of state for MGCG is p = —ap/ (1 + a) — $(z)p~ a /(l + a), 
where §(z) = -[po c (l + z) 3 ]( 1+tt )(l - n 0B ) a { a n aDM + n aDE [uj DE + a(l + u DE )](l + 
z j3uDi!(i+a) y Some cosmological quantities, such as the densities of different 
components of the universe il^ (i respectively denotes baryons, dark matter and dark 
energy) and the deceleration parameter q, are obtained. The present deceleration 
parameter go 7 the transition redshift zt and the redshift z eq , which describes the epoch 
when the densities in dark matter and dark energy are equal, are also calculated. To 
distinguish MGCG from others, we then apply the Statefinder diagnostic. Later on, the 
parameters (a and lode) of MGCG are constrained by combination of the sound speed 
Cj, the age of the universe to, the growth factor m and the bias parameter b. It yields 
a = — 3.07±|J| x 10~ 2 and uj de = — l-05±g;5i- Through the analysis of the growth of 
density perturbations for MGCG, it is found that the energy will transfer from dark 
matter to dark energy which reach equal at z eq ~ 0.48 and the density fluctuations 
start deviating from the linear behavior at z ~ 0.25 caused by the dominance of dark 
energy. 
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1. Introduction 

As standardized candles, type la supernovae suggest that the universe is undergoing 
accelerating expansion jl]. The same evidence has been shown by cosmic microwave 
background (CMB)|2j and large scale structure (LSS)[3] observations. All these data 
indicate the dominant component of the universe is relatively smooth and has a large 
negative pressure. All kinds of alternative models stand up for explaining this exotic 
phenomena within relativistic theory of gravity, which include the energy of the quantum 
vacuum (such as cosmo logical constant A[2J |5]), the existence of another new scalar 
field (quintessence [6] which is possibly related to the inflaton, k-essenceJ7], tachyon[8], 
phantom [9], and quintom[10]), and the influence of unseen additional spatial dimensions 
predicted by string theory [TTj. Instead of the presence of dark energy, the most 
interesting idea for solving this dark riddle is depend on a new aspect of gravity which 
is not accounted for general relativity[12]. Unfortunately, the current astronomical 
observations data still can not determine completely the nature of dark energy or decide 
its existence since this problem has puzzled us for a long time. 

As one of plausible dynamical models for dark energy, Chaplygin gas (CG) and the 
generalized Chaplygin gas (GCG)[l3], has created a lot of interest in recent times [T4"]. 
CG or GCG behaves like pressureless dust at early times and like a cosmological constant 
during very late times. The fact that the properties of CG and GCG interpolate 
between those of cold dark matter (CDM) and a A-term makes the models to provide 
a conceptual framework for a unified model of dark matter and dark energy. Besides, 
it is more worthy of note that the GCG model has been successfully confronted with 
various phenomenological tests such as high precision CMB data, supernovae data and 
gravitational lensingpj)]. On the other hand, since observations such as SNe lap] are still 
not sufficient to establish evidence of a dynamical equation of state for dark energy, we 
can naturally consider a fluid consisted of dark matter and quiessence dark energy with 
constant ujde which full of the whole universe. Thus, we propose a modified generalized 
Chaplygin gas (MGCG) as the unified dark matter and dark energy revisited on some 
previous works [T8| [T9| 120] . 

With future deeper and more intensive surveys of SNe la, such as Supernovae 
acceleration probe (SNAP)[29J, it provides a chance to distinguish different models of 
dark energy so that we will consider the so-called "geometrical " or Statefinder diagnostic 
since it can probe the expansion dynamics of the universe through higher derivatives of 
the scale factor. So, we must constrain the parameters (a and code) of the MGCG at 
first. The plat of this paper is as follows. In Sec. 2, we give corresponding fundamental 
cosmology equations about MGCG. Subsequently, in Sec. 3, we consider Statefinder 
diagnostic of the MGCG. We then constrain the parameters of the MGCG by means 
of the sound speed cj., the age of the universe to, the growth factor m and the bias 
parameter b in Sec. 4. Finally, conclusion and our results are outlined in Sec. 5. 
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2. Fundamental cosmology equations about MGCG 

2.1. Motivation 

An interesting unified model for two dark sectors is the Chaplygin gas (CG) model [16J 
introduced with p = —A/ p. Although this model has been very successful in explaining 
the SNe la data, it shows that CG model does not pass the tests connected with 
structure formation and observed strong oscillations of the matter power spectrum jTF]. 
This situation can be alleviated in the generalized Chaplygin gas (GCG) proposed with 
p = —A/p a . And the parameter a is rather severely constrained, i.e., < a < 0.2 at 
the 95% confidence level[l3]. Since the difference between ACDM and GCG models is 
so tiny and the equation of state of dark energy still cannot be determined exactly, a 
new generalized Chaplygin gas (NGCG)[20] proposed by considering that A in GCG is 
a function of redshift z. Following this idea, we develop GCG-like models from another 
view. As a purely kinetic k-essence model with a constant potential, the tachyon fields 
can be considered as CG model [IE]. Meanwhile the tachyon field can also act as a 
source of dark energy depending upon the form of the tachyon potential. Thus, the 
authors in Ref. [19] introduced an extended tachyon field (ETF) and then provided a 
modified Chaplygin gas (MCG)[19], in which the equation of state (EOS) has a more 
generalized form 

g 1 A m 

1 + a 1 + a p a 

where A and a are constants. As a kind of new attempt just like the extension from 
GCG to NGCG, we present a modified generalized Chaplygin gas (MGCG) through 
replacing A in Eq. ([I]) with a function d{z) of redshift z. The analysis and discussion 
of the matter power spectrum for MGCG will be investigated in our future work by 
comparing the model with observations. 

2.2. Deduction of some cosmological quantities and discussion 

Within the framework of Friedmann- Robertson- Walker (FRW) cosmology, we consider 
an exotic background fluid called modified generalized Chaplygin gas (MGCG) whose 
equation of state is as follows based on Ref.j 



a 1 "&(z) /nS 
p= p 2 

where a is a constant and d(z) is a function of redshift z. Since the MGCG is another 
new unified model with two dark sectors (namely, at early times the energy density 
behaves as dark matter: p oc (1 + z) 3 ; while at late times it behaves like a quiessence 
dark energy: p oc (l + z) 3(1+Wl3£; )), it could be supposed the energy density of the MGCG 
as follows 



P 



[«(1 + ^) 3(1+Q) + A(l + zf^+^ii+a)^ ; (3) 
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where the parameter of state for dark energy ojde is a constant and should be taken as 
any value in the range (—1.46, — 0.78) [21] . In Sec. 4, we will further constrain the lode 
in our model. Then the pressure of MGCG yields 

Ol\k{\ + Z) 3 ^ + A(l + 3)3(1+^X1+°)] + gfc) 
P ~ ~ (1 + a)[/«(l + Z )3(l+°) + A(l + 2 )3(l+^ B )(l+a)]a/(l+a) " U 

On the other hand, the whole pressure ptotai and energy density ptotai satisfy the 
conservation equation 

Ptotai + 3H(p total + p tota i) = 0. (5) 

According to conservation of baryons, the continuity equation is separated into two 
parts: 

p B + 3H{p B +p B ) = 0, (6) 
p + 3H(p + p) = 0, (7) 
where subscript "B" denotes baryons. Using the Eqs. (j3J), PJ and (j7|), we obtain 
= + a(l + u DE )]X(l + z )3(i+^^)d+«) 

+ a«(l + z) 3(1+Q) ). (8) 



Furthermore, the pressure of dark energy pde, the energy density of dark energy 
Pde and the one of dark matter p DM can be respectively written as 



Pde = P 

_ UJ DE \{1 + z)3(l+*g g )(l+tt) 

" + ^)3(l+«) + X(l + z )3(l+o; D B)(l+a)]«/(l+a) 

Pde = 

A(l + z )3(1+^db)(1+") 
" + z)3(l+a) _|_ X(l + z )m+u D E)(l+»)]a/(l+a) 

Pdm = P — Pde 

k(1 + ^) 3 ( 1+Q ) 



(9) 



(10) 



11) 



[AC(1 + ^) 3 ( 1+Q ) + A(l + 2 ) 3 ( 1+w -Di5)(l+«)]a/(l+") ' 

From Eqs. (PSJ and (03}, it gives that 

K + A = p^ +a = (poDM + P0DE) 1+a = Pl^ll - tt 0B ] 1+a , (12) 

= [k + A]«/(i+«)' (13) 

P° DE = [ K + A]^/M' (14) 
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where p oc , p Q , Podm and p 0DE are the present values of p totah p, p DM and p DE , 
respectively and the label "0" denotes today's evaluated quantities. Parameters k and 
A can then be written as 

Substituting Eq. (TI5]) to Eq. (jHJ), we have 

${ z ) = -Poc°[l " n B] a (l + z) 3(1+Q) {aftoDM 

+ £l 0DE [u DE + a{l + u DE )](l + Z f^+^Y (16) 
Then, the equation of state for the MGCG model reads as 




-1 1 2 3 4 5 6 

z 

Figure 1. The equation of state for the MGCG \( z ) as the function of redshift. 
Where n 0B = 0.0449, Qqdm = 0.222 and n 0DE = 0.734. 

The values of Qodm and Qode depend on the model used in the data precessing. 
For example, both the Wilknson Microware Anisotropy Probe (WMAP) |24j and two- 
degree field Galaxy Redshift Survey (2dFGRS) [25] used ACDM model of the universe. 
Besides, cosmological parameters are also constrained with other models. In this 
papar, the current density parameters used in the plots are Q b = 0.0449 ± 0.0028, 
Q QDM = 0.222 ± 0.026 and Q de = 0.734 ± 0.029 based on WMAP 7 data[24j. Later, 
we will give the constraints of a and to be i n section 4. From Eq.( lT7|) . we can plot the 
function x{ z ) as the function of redshift z ( see Fig. [1]). x( z ) is always negative from 
past to future. And we can see the effect of the parameter a and u E e on the x( z )- The 
evolution of x( z ) is mostly flat in the high redshift z > 1 and is very steep as redshift z 
becomes low. 



6 



1.005 



<x=-0.03, co DE =-1.05 
a=0.05, co DE =-1.00 
a=0.10, co DE =-0.98 




Figure 2. The Hubble parameter in units of Hacdm for the MGCG as the function 
of redshift. The priors floDM = 0.222, fl DE = 0.734 and Q b = 0.0449 have been 
used. 



When we consider a spatially flat FRW universe with the exotic background MGCG 
fluid and the baryon component, the Friedmann equation can be written as: 
8ttG 



Then, we obtain 



-ptotai- 



;is) 



E 2 (Z) EE 



IP 

m 



= (i - n 0B )Tf^[n 0DM (i + zf^ + n 0DE {\ + z f^n E )n +a)] ^ 
+ n 0B (i + z ) 3 . (19) 

The evolution of the Hubble parameter in units of -£/acdm is plotted in Fig. [2j During 
the cosmological evolution, the behavior of H is similar to -£?acdm at very early time, 
and when z — Y 0, H is approaching ACDM. 

By using the above equations, the densities of different components of the universe 
OflM, Qde and Qb can be respectively derived as 

(i-n 0B )^n 0DM (i + z) 3 



n 



DM 



n 



DE 



_ (l-^ 0B )T^^ 0Di ;(l + ^ 3 ( 1 +^) 

~ E\z)[n QDM {i + Z )-3« M (l+«) + fi 0D£ ]«/d+«) ' 

Oob(1 + ^) 3 
EHz) ' 



(20) 
(21) 
(22) 
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Figure 3. The densities of SIde, FIdm and £Ib as the functions of redshift. 

The relationships between fij (i respectively denotes baryons, dark matter and dark 
energy) and redshift z are shown in Fig. |3]with the prior of ^oda/ = 0.222, Qode = 0.734 
and Q b = 0.0449. Based on Fig. [3j the evolutions of baryons, dark matter and 
dark energy under MGCG model are consistent with what are recognized [211 E2]. An 
very important redshift described the epoch when the densities in dark matter and 
dark energy are equal is (1 + z eg )~ 3uJDE ^ 1+a ' ) = VIq DE /VLq DM - And z eq ~ 0.48 for 



a 



-0.03, u DE 



j eq ) 

-1.05, z eq ~ 



0.46 for a = 0.05, code 



-1.00, z, 



eq 



0.45 for 



a = 0.10, ojde — —0.98 respectively. From Fig. [3] we also know that the dominance of 
the dark energy leads to the acceleration expansion of our universe. The increasing 
density of dark energy is the reason why expansion of our universe transited from 
deceleration to acceleration. This transition will be described by our model in more 
detail with the deceleration parameter q. 
The deceleration parameter is 



q 



aH 2 



H 

IP 

■i + ^i^n 



x [a 



ODM 



2E 2 (z 

±{l + uj DE )n 0DE {l + zf" DE{1+a) 



OB 



(23) 



Thus we can give the following results: the present deceleration parameter q w —0.65 
according to z = and the transition redshift Zt ~ 0.70 according to q = when taking 
a = -0.03, u DE = -1.05, q « -0.60 and z T ps 0.65 for a = 0.05, u DE = -1.00, 
and go ~ —0.58 and zt ~ 0.60 for a = 0.10, ujde = —0.98, which are consistent with 
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observations. The deceleration parameter as the function of redshift is shown in Fig. HJ 
According to Fig. HJ the expansion of the universe is from slowing down (q > 0) in the 
past to speeding up (q < 0) at present time and in the future. And when q > is at 
the high redshift, q changes very slow with redshift z; when q < is at the recent time, 
and q changes faster and faster with redshift z. 




-1.2 1 1 1 1 1 1 1 1 

-1 1 2 3 4 5 6 

z 

Figure 4. The deceleration parameter q as the function of redshift. The priors 
il 0DM = 0.222, rt 0DE = 0.734 and O 0B = 0.0449 have been used. 

Now, we compare some cosmological characters between CG, GCG, NGCG, MCG 
and MGCG (see Tabled]). In order to distinguish these model, we use the same symbols 
as above. From Table [U we can see the EOSs of these models so that they have some 
different characters. First of all, CG and GCG can be treated as a unified fluid whose 
behaves like pressureless dust at early times and like a cosmological constant during 
very late times. For NGCG, it is considered as X-CDM, where X means quiessence, 
since the authors in Ref. [20] extended A of GCG model. And MCG is the mixture of a 
barotropic perfect fluid and a generalized Chaplygin gas. At the same time, considering 
A of Eq. ([[]) changed with redshift z, we proposed MGCG. 

Among them, GCG is a extended model of CG. When a = 1, GCG model becomes 
CG. And NGCG is a extended model of GCG. When A(z) = const., NGCG reduces to 
GCG. MGCG is a extended version of MCG. When ${z) = const., MGCG yields MCG. 
Meanwhile, both GCG and MCG can be conceived as ETF, in which GCG corresponds 
to the simplest model driven by a constant potential [19J. It is worth of noted that CG, 
GCG, NGCG, MCG and MGCG are all interact models. Namely, there are energy 
exchange between dark matter and dark energy. Basically, the behaviour of the energy 
exchange can be presented as Pdm/pde- We can see the transfer direction of the energy 
flow for CG is from dark matter to dark energy. However, the transfer direction of the 
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Table 1. Compare CG, GCG, NGCG, MCG and MGCG 



Models EOS A Parameter(s) Pdm/pde 



CG 


P = 


_A 

P 




Const. 


A 


~(l + z] 


i 3 


GCG 


P = 


A 

re 




Const. 


A and a 


~(l + z) 


|3(l+o) 


NGCG 


P = 


A(z) 




Function of z 


ude and a 


~(l+z) 


j-3wDE(l+a) 


MCG 


P = 


-ThP- 


1 A 
1+a p a 


Const. 


A and a 




|3(l+a) 


MGCG 


P = 




1 ti(z) 
1+a p a 


Function of z 


ude and a 




j-3wDB(l+a) 



energy flow for the others only depend on the parameter a if a ^ 0. When a > 0, the 
transfer direction of the energy flow is from dark matter to dark energy. When a < 0, the 
direction of transfer is inverse. Thus, constraints of parameters such as a and ujde from 
observational data and theoretical analysis will play an important role in understanding 
the physical nature of these models. For GCG, authors in Ref. [27] obtain a < 10 -5 
or a > 3 (the sound speed is larger than the speed of light) by using gauge-invariant 
analysis of perturbation. For NGCG, constraints of parameters are 1 + a = 1.06io!i6 and 
code = — 0.98^o;2o by observational constraints from SNe la, CMB and LSS data[20j. 
In section 4, we will constrain a and iode in our model by the maximum likelihood 
estimation. 



3. Statefinder diagnostic 

Dark energy models such as the cosmological constant [U [5], quintessence [6], K- 
essence[7J, Chaplygin gas [13] and quintompIJ] etc. have properties which can be 
model-dependent. In order to distinguish the very distinct and competing cosmological 
scenarios involving dark energy, a sensitive and robust diagnostic of dark energy is 
needed. So a new diagnostic of dark energy called "Statefinder" diagnostic has been 
constructed by Sahni et al. [23] who were using both the second and the third derivatives 
of the scale factor a. 

The statefinder pair {r, s}, in addition to the oldest and most well-known geometric 
variables H and q, defines two new cosmological parameters 

' s 7Jp- < 24) 

r — 1 

s = — — . (25) 

3(9-1/2) 1 J 

An important property of the Statefinder is that spatially flat ACDM corresponds to 
the fixed point 

{t-,s}|acdm= {1,0}. (26) 

Clearly an important requirement of any diagnostic is that it permits us to tell difference 
between a given dark energy model and the simplest of all models - the cosmological 
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constant just as demonstrated in[23j-[28j. By using the r(s) evolution diagram, the 
discrimination between a given dark energy model and the ACDM scenario can be clearly 
identified. And it is more worth of noted that the Statefinder diagnostic combined with 
future SNAP [29] observations may possibly be used to distinguish between different 
dark energy models. 

Based on the Eqs. ( ]2lj) and fl25l) . we obtain the Statefinder parameters for the 
MGCG as follows 



1 + 3 + = 1 + 3T(s) + 6(z) 



6H/H 2 + 2H/H' 3 
9 + 6H/H 2 



6T(z) + 2Q(z) 
9 + 6T(» 



(27) 
(28) 



where 



T(z) = 



H 

IP 

3 

2 



-udeQ 



ODE 



(1 + z f^ E (l+a)^ _ Q 0B )«/(1+«)|(1 _ Q oB )a/(l+a) 



X 



ODM 



ODE 



'l + Z 



OB 



n 



ODM 



ODE 



[l + Z 



,3u) DE (l+a) 



itt/(l+a)l -1 



(29) 



and 



2 2 

+ (i + a )^ E n ofl£ (i + ^ (1+a) 



ODE 



'l + Z 



>3u) DE (l+a) 



-1 



ODM 



n 0B ) a/{1+a) 



X 



n, 



ODM 



OB 



n 



ODM 



ODE 



+ + z 



ODE 



'l + Z 



3u) DE (l+a) 



a/(l+a)i -1 



(30) 



Then, the evolution trajectories of our model in the plane of the Statefinder parameters 
can also be plotted (see Fig. [5] and Fig. [6]). 

Since different dark energy models have distinct evolution trajectories in the plane of 
the Statefinder parameters, we can distinguish between various dark energy models and 
the MGCG. Just as the literatures |23[ 128] show: the ACDM scenario corresponds to the 
fixed point {r = 1, s — 0}; the standard cold dark matter (SCDM) scenario corresponds 
to the point {r = l,s = 1}; the Statefinder parameters plane of quiessence are some 
vertical segments where r decreases monotonically from 1 to 1 + %ojde(1 + <^de) and s 
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Figure 6. The Statefmdcr r(s) evolution diagram when we fix oj de and vary a. Dots 
locate today's values and arrows denote the direction of evolution for the MGCG. 
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remains constant at 1 + ude\ the quintessence tracker models have typical trajectories 
similar to arcs of an upward parabola lying in the regions {r < 1, s > 0}. 

Now, we can see the Statefmder parameters plane of the MGCG from Fig. [5] and 
Fig. E In Fig. El we fix a = 0.5 and vary u DE as -1.05, -1.00, -0.95 and -0.90 
respectively. Where ude = —1-00 exhibits a complete downward parabola, while the 
others correspond to some broken parabolas. And from this diagram, we can see the 
Statefmder trajectory begins with r — 1. Besides, the cases of ujde < — 1 never arrive at 
the fixed point {1, 0} and when ude > —1 the trajectories must pass through this point. 
Finally, we can see that the less ude, the smaller today's value s and the larger today's 
value r . However, in Fig. |6j we fix ude = —0.85 and vary a as 3.1, 2.7, 2.3, 1.9, 1.51, 
1.1 and 0.7 respectively. It is interesting to see that the trajectories can pass through 
the fixed point {1, 0}. And when fixed code, the less a, the smaller today's value sq and 
the larger today's value r . We notice the present Statefmder points as well as the fixed 
point {1,0} locate on a straight line which means the relationship between r and Sq is 
linear when we fix ude- That is because the following relationships 
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r o = 1 + -j^> DECODE 



-,,/-,, n OiLJ DECODE 

1 + (1 + a)u DE 



n 



0B 



so that 



s = 1 + (1 + a)u DE , (31) 

l — \l 0B 

9 

r = 1 + -ude^odeSq- (32) 

Through above, the separation between distinct families of dark energy models is 
very remarkable when we analyze evolutionary trajectories using the Statefmder pair 
{r, s}. Forthcoming space-based missions such as SNAP are expected to greatly increase 
and improve the current Type la supernovae inventory, and maybe can identify which 
kind of dark energy model. 



4. Constraints on MGCG 

4-1. The sound speed 

Since our model can be considered as a new scenario for unified dark matter and dark 
energy, we have to study its density perturbations and the structure formation, we also 
have to investigate its adiabatic sound speed. So, in this section, we consider some 
effective parameters for the MGCG such as a and ujde by the maximum likelihood 
estimation. If the MGCG is considered as a perfect fluid satisfying equation (J2]), then 
the MGCG component will cluster gravitationally with the adiabatic sound speed given 
by 
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Figure 8. The region inside the curve is the domain of parameters a and ujde which 
is constraints from the sound speed. 
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a _ a tt 0DE {u DE + g(l + u DE )](l + zgWi+j + gfl ogM 
1 + a 1 + a n 0DE (l + z) 3uJ o E (i+ a ) +n 0D M 

(1 + a; gE )Oo DE [a; D£ + a(l + a; gJS )](l + gfWi+j + a fl ogM 

From the above equation, we obtain a / -1. However, when a < —1, we can 
derive (? s \ z -+oo — <^de < 0, and when c 2 is negative, the MGCG fluid is unstable. When 
a > —1, we have c 2 \ z ^oo ~ 0. Thus, the parameter a must be greater than —1 because 
the limit of the sound speed was very small at early time and can not be larger than 1 
even though in the future if our model can be as the unified dark matter-dark energy 
From Fig. \J\ we can see the slope of the curves for c 2 depends on a and uj ee . Besides, 
we clearly see that deep in the matter era the behavior of the MGCG closely resembles 
that of CDM. 

To constrain a and to be further, combined with the data of SNe la, CMB and 
2dFGRS, it imposes a range on ube- —1.46 < code < — 0.78 [21J . Through the limit of 
the sound speed < c 2 < 1, when we use the Eq. ( 155]) . we can find the domain of a - 
ujde plane (see Fig. [8]). Then, we derive —0.3 < a < 4.7 and —1.19 < ujde < —0.78. 
From Fig. [8] we can see the boundary of domain is constraints from the sound speed. 



4-2. The age of the universe 

By using a distance-independent method, Jimenez et al.[32] determined the age of 
globular clusters in the Milky Way as to = 13.5 ± 2Gyr. By using the white dwarfs 
cooling sequence method, Richer et al. [33] and Hansen et al.[3l] constrained the age 
of the globular cluster M4 to be to — 12.7 ± 0.7Gyr. Then, the age of the universe 
needs to satisfy the lower bound: to > 11 — 12Gyr. Assuming ACDM, WMAP7 data 
give the age of the universe t = 13.73 ±0.13Gyr[2l]. Here we adopt model-independent 
astronomical observations of globular clusters as a criterion. Making use of the definition 
for the Hubble parameter H = a/a and the relationship between the scale factor and 
the redshift 1 + z = 1/a, then we have 

da 1 , 

*^=W 2 ' (34) 

Integrated the above equation, we obtain 

tn= dt= — 

Jo Jo H(l + z) 

1 r°° dz 



Ho J0 (1 + 2)5/2^(1 _ n QB )TTZ{n QDM + Q 0DE (l + 2 )3aWl+a)]Tfe + Q 0B 

In Fig. EH with using the age data of the universe, we plot the isochrones in a - u^e 
plane under the boundary of domain constrained by c 2 . Then, we can see that the 
parameters a and uode are constrained by the age of the universe further. Since the 
luminosity distance for all models is given by the simple expression for a spatially flat 
universe 

d /Driz^T- 1 



Hz) 



•v 

dz V 1 + z 



(36) 
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using the Eqs. (I35]) - (15S|) . we have 



1 + Z 



.— cte. (37) 
o Jo (l + z) 3 v ; 

From above equation, we can see that the luminosity distance is directly related to the 
age of the universe. 
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Figure 9. The dashed-dot curve is the domain of parameters a and lude which is 
constraints from the sound speed. The dot lines is the isochrones for parameter to in 
the a-uiDE plane. 



4-3. The linear density perturbations 

In this subsection, we will study the growth of density perturbations for the MGCG fluid 
in the linear regime on subhorizon scales. We will follow the analysis method in Ref . |13] . 
Firstly, we will introduce some quantities used in this part. The relationship between the 
comoving coordinate (x, t) and the physical coordinate (r, t) is x — f/a. The velocity 
in physical coordinate is u = r. Furthermore, we can expand u = ax + v where v is 
the first-order perturbation to the Hubble flow ax. pDM(x,t) = poDM(t)(l + o~DAi(x,t)) 
where Podm(^) is the background density and SjjM{x,t) is the first-order perturbation 
of Podm- 4>{x, t) is the first-order perturbation of Newtonian potential $, which satisfies 
Poisson's equation V 2 $ = 4nGa 2 p£,M in the comoving coordinate. 

Before discussing the linear density perturbations, we will research an explicit 
interaction between dark matter and dark energy in our model and depict this interaction 
through an energy exchange term T which will be used in the linear density perturbed 
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equations. From Eqs. ( llOp - (TTTj) . we can obtain the scaling behavior of the energy 
densities as follows 

P DM _ 1m _j_ z \-3u DE (l+a) _ MqDM ^ _|_ ^-3 WI)E (1+ Q ) 

Ade A Qqde 
where there is an explicit interaction between dark matter and dark energy. This can 
be seen from the energy conservation equation more clearly, which can be written as 

Pdm + 3Hp DM = [pde + 3H(p DE + p DE )). (39) 

Furthermore, 

Pdm + 3Hp DM = T, (40) 

p DE + 3H(p DE + p DE ) = -T, (41) 

where T is the source term in the continuity equation due to the energy transferring 
between dark matter and quiessence dark energy in our model. From the above 
equations, we can respectively derive the effective equations of state for dark matter 
0Jdm and dark energy u e Jl in the MGCG scenario 

r 



eff 
DM 



3Hp D M 
auD E n Q DE{l + zf^ l+ ^ 

'n 0DM + n DE(l + z) 3 ^ 1+a) ' 



(42) 



eff , 
U DE — U DE + 



3Hp 



DE 



= " DE(1 + a) -?wT?wTTIF^- (43) 

Now, we focus on the mass density perturbations in the linear approximation. 
In this paper, we apply the comoving coordinate (x, t) for a pressureless fluid with 
background density Pqdm and density contrast 5dm, with a source term Y due to the 
energy transfer from dark matter to quiessence dark energy. Based on Ref. [30], in 
the first-order approximation, the perturbed equations are written in the comoving 
coordinate as follows 

-j- + -v = ^, 44 

at a a 

„ 95dm r S D M 

V • v — — a— ar , (45) 

<ft PODM 

i-V 2 = AttGpodmSdm, (46) 
a z 

where 5dm, v and are respectively the first-order perturbations of Podm, <ix and the 
Newtonian potential $. Since the matter is so cold that the term proportional to the 
speed of sound c 2 s in Eq. ( 144"]) can be omitted. It leads to disappearance of the term 

2 

V5da/ in the Eq. (1441 [3Tj . Assuming that both the density contrast 5dm and peculiar 




Figure 10. 5dm as the function of scale factor a. 
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Figure 12. The bias b(a) as a function of scale factor a. 




Figure 13. The contours for growth factor m in the a - lode plane. 
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Figure 14. The contours for bias b in the ujde - a plane. 




Figure 15. The contours in solid lines come from four limits of c 2 a , t 0} m and b in the 
lode - a plane. The confidence level contours of 68.3%, 95.4% and 99.73%. The la fit 
values for our model parameters are a = — 3.07+4gg x 10~ 2 and lu de = — 1.05l ; n . 
And the plus sign in this figure denotes the best fits. 
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velocity v are small, i.e., that 8dm ^ 1 an d v <C u, where u is the velocity in physical 
coordinate. Taking the divergence of Eq. ( jSJ) and substituting Eqs. f ]4"5]) and fH6|) . we 
finally obtain 



EH" DM +(h+ 2H 2 + H )s 

V Pdm J 



DM 



r / r x/n 

2# iff 

Pdm \Pdm 



AttGpdm 

dm = 0, (47) 



where ' = d/d In a. Since the term ^V8dm is not included in the Euler equation, it 
leads to disappearance of the term V 2 8dm in Eq. (14Tj) . And we can easily see that there 
is no scale dependent term to drive oscillations or blowup in the power spectrum. We 
can solve Eq. (H7|) numerically, when the initial conditions are chosen as a = 10" 3 and 
8dm — 10 -3 , we plot the linear density perturbation for dark matter 8dm as a function 
of a (see Fig. [TOl) . In Fig. [TUl we can see the perturbation starts departing from the 
linear behavior around a ~ 0.80 i.e. z ~ 0.25 which is similar to the epoch when z < z eq . 
Through the energy transferring from dark matter to dark energy, we can see that the 
dominance of dark energy is related to the time when the density fluctuations start 
deviating from the linear behavior. 

When we consider the period after decoupling, the baryons are no longer coupled 
to photons so that the baryons are also a pressureless fluid like the dark matter. Then, 
we derive the baryon perturbation as follows 

R 2 8" B + (H + 2H 2 )8' B - AttGpdmSdm = 0. (48) 

Then we investigate the growth factor m{y) and the bias b(a) where m(y) = D'(y)/ D(y), 
D(y) is the linear growth function, y = ln(a) and b = 8b /8dm- One can see from Fig. [TT] 
and Fig. [12]that between the present and z ~ 5, the growth factor m(y) and the bias b(a) 
are quite sensitive to the values of a and ude- Subsequently, we consider the constraints 
from m(y) and b(a) on the a - ude plane. In Fig. [13] and Fig. [TH we have shown the 
contours for the growth factor m(y) and the bias b(a) in the a - iode plane under the 
boundary of domain constrained by c 2 . On the other hand, the growth factor and the 
bias parameter at z ~ 0.15 have been determined using the 2dFGRS. The redshift space 
distortion parameter /3 = 0.49 ±0.09, and the bias b = 1.04 ± 0.11 [35]. For /3 = m/b, we 
can subsequently determine the constraint on the growth factor m as m = 0.51 ± 0.11. 
Noted that the fitted values of b and (3 are sensitive to the underlying model and the 
observational constraints on them are obtained under ACDM when converting redshift 
to distances for the power spectra. The right way is to conduct data processing for b 
and (3 by using 2dFGRS data under our model. Then we can constrain a and ujde by 
using the fitted values. In this paper we use the above fitted data temporarily and our 
future work will focus on this tough work. 

In Fig. [T5J we consider four restrictions together to give the best estimation of a 
and ude- Firstly, we plot contours constraint imposed by to (Fig. [9]), m (Fig. [13]) and 
b (Fig. H3]) under the boundary of domain constrained by the sound speed. Each set 
leads to a reduced x 2 : x 2 ( a > <^de) =(Theory values-Observation values) 2 / (Observational 
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error) 2 . To combine the constraints on a and U) EE coming from the age of the universe to, 
growth factor m and the bias b within the domain given by the limits of the sound speed 
c 2 , we have added their individual x 2 as if they were part of a total experiment with 
uncorrelated Gaussian errors: x 2 totai( a ^DE) = X 2 (^,^de) + Xm(oi,u DE ) + xl(a,u DE ). 
Therefore, we plot the contour level Ax 2 ota i(a, ojde) = 2.3, Ax 2 ota i(a,ujDE) = 6.17 
and Ax 2 tota i(a,u DE ) = 11.8, where Ax 2 otal (a,u DE ) = x 2 (a,u DE ) - (x 2 (a, u; DE )) TOi „, 
defines, respectively, for two degrees of freedom the 68.3%, 95.4% and 99.73% confidence 
levels represented in Fig. [151 The la fit values are a = — S.C^l^gg x 1CT 2 and 
uj ee = — 1.05jlo!ii- Noted that the presence of the cutoff contours' curves in Fig. [T5l is 
due to the restricted area boundary which comes from the limits of c 2 . 



5. Conclusion 



A new model named as the modified generalized Chaplygin gas (MGCG), which is 
a further generalization of the generalized Chaplygin gas, has been proposed in the 
present paper. As a version of the unified dark matter and dark energy, this MGCG 
fluid is consisted of dark matter and quiessence dark energy with constant u EE . Firstly, 
fundamental cosmology equations for the MGCG have been described. For tests with 
future deep observations, we then consider the Statefinder diagnostic since it can probe 
the expansion dynamics of the universe through higher derivatives of the scale factor. 
Furthermore, we investigate the evolution of density perturbations and the structure 
formation in our model. Then, with applied the age of the universe data t , growth 
factor m and the bias b within the domain given by the limits of c 2 , the parameters a and 
u EE are constrained: a = — 3.07lf;gg x 10~ 2 and uj de = — 1.05io!ii- It has been shown 
that there are a little difference between our model and others. To determine whether 
our model could be a candidate model of dark energy, the analysis and discussion of 
the matter power spectrum for MGCG must be done. This is a subject of our future 
research. 
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